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Control of magnetism in singlet-triplet superconducting heterostructures 
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We analyze the magnetization at the interface between singlet and triplet superconductors and 
show that its direction and dependence on the phase difference across the junction are strongly tied 
to the structure of the triplet order parameter as well as to the pairing interactions. We consider 
equal spin helical, opposite spin chiral, and mixed symmetry pairing on the triplet side and show that 
the magnetization vanishes at 0 = 0 only in the first case, follows approximately a cos 0 behavior 
for the second, and shows higher harmonics for the last configuration. We trace the origin of the 
magnetization to the magnetic structure of the Andreev bound states near the interface, and provide 
a symmetry-based explanation of the results. Our findings can be used to control the magnetization 
in superconducting heterostructures and to test symmetries of spin-triplet superconductors. 


Introduction. Attempts to combine dissipationless 
transport in superconductors with the control of spin cur¬ 
rents has driven many recent studies of superconducting 
heterostructures [T]. In the vast majority of superconduc¬ 
tors the conduction electrons pair in a spin-singlet state, 
freezing out the spin degrees of freedom at low temper¬ 
ature. However, triplet correlations near interfaces with 
magnetic materials produce equal-spin Cooper pairs that 
can propagate through a magnet leading, for example, to 
the long-range proximity effect las]. 

The alternative pathway to spin control by supercon¬ 
ductivity is via utilizing compounds that support spin- 
triplet pairing. The number of known triplet super¬ 
conductors (TSCs) has been growing steadily, and now 
includes UPts [5], ferromagnetic superconductors such 
as UGe2, URhGe, Ulr and UCoGe [SHH], quasi one¬ 
dimensional organic system (TMTSF)2X (X=C104 and 
PFe) [21110]. Singlet and triplet states are mixed if the 
material lacks inversion symmetry mills, and among 
the non-centrosymmetric superconductors Li2Pt3B has 
the clearest indication of a significant triplet compo¬ 
nent [13 ns. The strongest evidence for triplet supercon¬ 
ductivity has emerged for Sr2Ru04 mum. Existence of 
very pure single crystals with the perovskite structure 
made this material a testbed for studying heterostruc¬ 
tures based on triplet superconductivity miiiH]. 

From energetic considerations most of the non¬ 
magnetic superconducting compounds that support 
triplet pairing should have a unitary order parame¬ 
ter [19], so that the Cooper pairs do not have a net av¬ 
erage spin-derived magnetic moment m- Hence, triplet 
superconductors cannot by default be assumed to sup¬ 
port dissipationless spin transport. However, in con¬ 
junction with superconducting orders of different sym¬ 
metry, nontrivial spin aspects of triplet superconduc¬ 
tivity appear. Andreev bound states in singlet-triplet 
Josephson heterostructures were analyzed [2T1124] . and 
spin-accumulation was found when a phase difference 
is established across the junction mHaZ]. In parallel, 
it was shown that admixture of the sub dominant order 
may lead to spin accumulation, spin and charge currents 
near a boundary of a chiral triplet superconductor [28] 


or anomalous flux response in mesoscopic loops [29] . 

In this Letter we show that the magnetism in singlet- 
triplet superconducting heterostructures is fundamen¬ 
tally linked to the nature of the triplet pairing. We con¬ 
sider a microscopic model for a high-transparency inter¬ 
face between a singlet and a unitary triplet (or mixed par¬ 
ity, see below) superconductor, and self-consistently solve 
the corresponding Bogoliubov-de Gennes (BdG) equa¬ 
tions to obtain the energy spectrum. In all cases we 
find spin-splitting of the Andreev bound states (ABS) 
near the interface, leading to a magnetization parallel to 
the spin-triplet d-vector. However, the variation of the 
magnetization M with the phase difference (j) across the 
junction is non-trivial, and depends on the nature of the 
pairing interaction on the triplet side. Our main results 
are summarized in Fig. I, where we consider three dis¬ 
tinct cases. If there is no singlet component of the pairing 
interaction on the triplet side, M oc sin(/) only appears 
under a finite phase difference, in agreement with Refs.l^ 
and |26l This is realized, for example, for equal spin 
pairing states. Real space interactions leading to triplet 
5'^ = 0 pairing often promote a (subdominant) singlet 
pairing, so that the singlet amplitude persists into TSC 
near the interface. The magnetization resulting from this 
mixed symmetry, M oc cos0, already occurs at 0 = 0. 
Finally, motivated by the results of Ref. |28l we consider 
a mixed parity superconductor in contact with a singlet 
counterpart, and show that the magnetization has a com¬ 
plex dependence on (f. 

These results show that the phase difference can be 
used to control interface magnetization. In addition, the 
interface magnetization probes the spin structure of the 
pairing interactions. This finding is especially relevant 
for Sr2Ru04. There is no consensus on the exact form of 
the triplet pairing in this system, and NMR Knight shift 
measurements 1301 , absence of edge currents m, and the 
observation of half-quantum vortices [32] put in doubt 
strong pinning of the triplet order spin vector, d{k), to 
the crystalline c-axis. The chiral, d{k) = z{kx-\-iky), and 
the helical, d{k) = xkyFyk^^ states compete closely, [33] 
and distinguishing between these two possibilities with 
similar bulk gap is an important application of our re- 
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FIG. 1. Schematics for the superconducting heterostructures 
between spin-singlet (SCI) and triplet-active (SC2) super¬ 
conductors. Panel a): helical equal spin pairing state with 
no pairing interaction in the singlet channel. Panel b): chi¬ 
ral pairing with finite interaction in the singlet channel (see 
discussion in text). Panel c): mixed-parity order parameter 
near the interface. Vsi (Vn) are the couplings in the sin¬ 
glet (triplet) channel for sides z = 1,2. We show the profile 
of the order parameters (full line), with dot-dashed line for 
the proximity-induced order. We also indicate the qualitative 
dependence of the interface magnetization on the phase 0. 


suits. 

Model and formalism. We consider a two-dimensional 
lateral heterostructure made of two superconductors, 
SCI and SC2, having spin-singlet (SCI) and spin-triplet 
or mixed-parity (SC2) pairing. The x and y planar di¬ 
rections are perpendicular and parallel to the SC1/SC2 
interface, respectively. The system is uniform along the 
y axis, so that the translational symmetry is broken only 
in the x direction. The Hamiltonian is then defined on a 
square lattice of size Lx L (the lattice constant is unity), 
with periodic boundary conditions along y, 

H = Hi LH 2 L Hi2 (1) 

with {m = 1,2) being the Hamiltonians of each su¬ 
perconducting side 

Hm — ^ ^ (^i cr^j h.C.) fljn ^ ^ ^icr (2) 

~ Kn ^icT^jcr' ~ ’ 

(iJ> i 

and the interface term, 

Hi2= X] + h.C.) . (3) 

<5=±1 


Here the lattice sites are labelled by i = {ix^iy}^ with 
ix and iy integers between —1//2 and I//2, (i,j) denote 
nearest-neighbor sites, and y is the chemical potential. 
Labels 0 = {Og?/} and ±1 = {±l,Zy} denote the sites 
at the interface and their nearest-neighbors, respectively. 
The attractive interaction — (Kn > 0) can be chosen 
to be effective in the Sz = 1,0,—1 projections for the 
the TSC and/or in the S = 0 singlet sector. The local 
attractive term —Um {U > ^) only promotes spin-singlet 
pairing. For simplicity we take ti = t 2 = tj_ = t, and use 
the hopping parameter t as a unit of energy. Below, in 
the description of the results we discuss the qualitative 
consequences of relaxing this assumption. 

To investigate the model of Eq. 0 we decouple the 
interaction term in the Hartree-Fock approximation by 
introducing the pairing amplitude on a bond, = 

(cicrCjcr'), and on-site Aq = so that 




. (j(j' 12 \ 


These expressions yield the spin singlet {S) and triplet 
(T) components in the Sz = 0 sector, ± 

AJj^)/ 2, and the triplet pairing in the Sz = {1,-1} sec¬ 
tors, /MfS = A^^. They define in turn the supercon¬ 
ducting pair amplitudes with s- or p-wave symmetry, i.e. 
A,(i) = (A5+, + A5_, + A^,^, + A^,_,)/4, = 

- ^A-i(y))/2 and which 

are then determined self-consistently [34] . It is important 
to note that, while the terms generate equal spin 

triplet pairing, the coupling generically promotes 
both triplet 5'^ = 0 and non-local singlet (extended 5- 
wave) pairing. Similar decoupling of the U term produces 
local singlet pairing only. 

By suitable choices of the interactions and 1/ we 
can therefore model interfaces between superconductors 
with different symmetries of the order parameter. We 
take the most conventional SCI with 1/ 0, = 0, 

and compare the results for three different choices of SC2. 
Case a) is equal spin triplet pairing, ^ 0; = 

U = 0. In the bulk this corresponds to the helical state, 
A(A;) = i{d{k) • cr)cr^ with d{k) = xky + ykx^ where 
A{k) = f drijAij exp(i/c -rij) and Aq is the order 

parameter in spin space defined above. Case b) is a TSC 
with 5'^ = 0 and a subdominant extended s-wave pairing, 
yn=yU ^Q.yU ^ 0. Finally, case c) is a TSC 
with S^ = 0, 7 ^ 0, but possible local 5-wave pairing, 

U ^0, while = 0. 

For the last two situations the most favorable pairing 
state due to the depends on the electron density, 
n, and the chiral kx + iky order is stabilized in the re¬ 
gion between low doping y cx 1.2 and high (low) den¬ 
sity {\y\ cx 2.25) [35]. Hence we choose \y\ cx 1.8. All 
the numerical results below have been obtained for non¬ 
vanishing components of the pairing interaction V = 2.5, 
and U = 2.5, and system size L = 120. Greater values 
of L and modification of the couplings leave the results 
qualitatively unchanged. To investigate the effects of the 
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FIG. 2. Order parameter and magnetization for heterostruc¬ 
tures between a singlet (left side of the junction) and different 
triplet systems. Panels (a) are for helical equal spin pairing, 
(b) for chiral = Opairing, and (c) for the mixed symme¬ 
try case. Left column shows the evolution of the singlet and 
triplet pairing amplitudes, while the right column shows the 
net magnetization along the direction of the d-vector compo¬ 
nent for which the surface is pairbreaking as a function of the 
phase difference, 0. 


phase difference between the two superconductors, as is 
commonly done for the Josephson junctions, we trans¬ 
form the pairing wave-function in the SCI and SC 2 by 
the phase factors exp[—i 0 / 2 ] and exp[i 0 / 2 ], respectively. 
Self-consistent solution of the BdG equations gives the 
spin-resolved energy spectrum of the system, including 
both bulk and the surface (Andreev) states. The spec¬ 
trum determines the surface magnetization, and we com¬ 
pute M by summing the spin expectation values over all 
the occupied states at zero temperature. 

Results and discussion. The interface breaks the trans¬ 
lational symmetry in the x direction, and is therefore 
pair-breaking for the component of the triplet pair¬ 
ing. The differences between the three investigated cases 
stem from the consequences of this suppression for the in¬ 
terface magnetization. Case a) corresponds to the usual 
proximity coupling with step-like change of the symmetry 
of the pairing interaction across the junction, where the 
reduction of the amplitude of the Px component simply 
changes the effective coupling between the two sides, see 
Figs. 2 (a). The magnetism of the Andreev states local¬ 
ized at the interface is due to the splitting of the energy 
levels for opposite spins. This splitting appears because 


the phase shift between the singlet and the triplet en¬ 
ters the continuity condition for the wave function at the 
boundary with opposite signs for the two spin orienta¬ 
tions. Hence the magnetization vanishes at 0 = 0, and is 
nearly constant away from this point, see Fig. 2(a2). The 
situation is analogous to that discussed in Ref. |26l where 
a step-like change was found without the self-consistency 
on the order parameters, and assuming a delta-function 
barrier at the interface (BTK approximation). 

The origin of this dependence is easy to understand 
from the Ginzburg-Landau (GL) expansion of the free 
energy density, with the relevant terms / = aw? + 
ibm{'ilj^r]x — Here a > 0 is the inverse suscepti¬ 

bility, (rjx) are the amplitudes of the singlet (triplet 
kx) order, and m is the component of the magnetization 
along the d-vector coupled to kx^ in our case y. Near the 
interface where both -0 = and = \r]x\e'^^^‘^ co¬ 

exist, linear coupling to the magnetization ensures that 
the mimimum of the free energy is at m oc sin (j). The GL 
analysis is valid near Tc where indeed a sin 0 -like shape 
develops from the step-like behavior of the same symme¬ 
try. [26] The spatial profile of the magnetization, shown 
in Fig. 3(a), confirms that it is constrained to the region 
of the order of the coherence length around the junction’s 
interface. 

In general the magnetization linearly couples to the 
component of the d-vector for which the surface is pair 
breaking, and hence will be along the £ axis for the re¬ 
maining two situations. The phase dependence also be¬ 
comes more complex if there is pairing in the subdom¬ 
inant s-wave channel on the triplet side, as in cases b) 
and c). Recall that we have a chiral triplet order stabi¬ 
lized in the bulk. Reduction in the triplet kx component 
near the interface generically implies 7 ^ — ATj_^, 

and allows for admixture of singlet pairing. The dif¬ 
ference from the well-known case of sub dominant pair¬ 
ing near the surfaces in d-wave superconductors m is 
that it is the same coupling constant that promotes 
coupling in both channels, and therefore the admixture 
is at least parametrically stronger than in the case of 
sub dominant coupling. Indeed, it was found that even 
small variations in the surface barrier at the boundary 
of a triplet superconductor generated a substantial ad¬ 
mixture of a singlet component and emergence of mixed 
parity superconductivity near the edge [28]. Assuming 
that it is the subdominant extended 5-wave that is cou¬ 
pled to the triplet component near the interface, the lin¬ 
ear coupling terms in the GL expansion take the form 
fm = b'm[ip{xydxTj{x) + 'ip{x)dxTj{xy], and promote 
m oc COS 0 phase dependence of the magnetization (a 
similar term was found in the microscopic theory [36]). 
This is illustrated in Figs. 2(bi) and 2(ci), where the 
ky component of the chiral order is suppressed concomi¬ 
tantly with the kx component (in contrast to Fig. 2(ai)), 
and the magnetization already exists without the phase 
difference across the junction. 

In principle both sin0 and cos0 contributions should 
be present in such a junction, leading to a non-trivial 
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phase dependence. However, we believe that in case b) 
the competition with the subdominant pairing on the 
SC2 side causes the triplet components to decay more 
rapidly towards the interface and save energy via singlet 
pairing (compare Fig. 2(ai) and Fig. 2(ci)), enhancing 
the latter term, and exhibiting dominant cos (p behavior. 
Indirectly this is also supported by the magnetization 
that is more sharply localized near the interface, com¬ 
pare Fig. 3(a) and Fig. 3(b). 

Finally, inclusion of the on-site s-wave pairing for case 
c) dramatically extends the range of coexistence of triplet 
and singlet pairing in SC2, see Fig. 2(ci). This implies a 
more pronounced competition between the non-gradient 
and gradient couplings of the superconducting orders to 
the magnetization, but there is an additional complex¬ 
ity because the “local” s-wave component suppresses the 
ky component of the pairing faster than the compo¬ 
nent, Fig. 2(ci). The induced 5-wave order on the SC2 
side is phase-locked to the triplet component, yielding 
a finite magnetization that only weakly depends on the 
phase for p G (0,7r/2), but changes sign as the phase 
difference approaches tt. Presumably this occurs because 
of additional states appearing due to the singlet com¬ 
ponent changing sign across the junction. This view is 
supported by considering the spatial profile of the mag¬ 
netization, Fig. 3(c), which shows a two-peak structure 
for small 0, including the contribution away from the in¬ 
terface, but only exhibits a single peak at the interface 
for 0 = TT. The resulting phase dependence of the mag¬ 
netization contains many Fourier components. 

Analysis of the energy spectrum. Fig. 3(d), supports 
these conclusions. If the structure of the order parameter 
near interface plays the major role in the emergence of 
the magnetization, nearly grazing trajectories must have 
a significant contribution. Indeed, the spin-split branches 
of Andreev states crossing the Fermi surface are moved 
to large values of /c^/tt ^ 0.5. In addition, there is a sec¬ 
ondary branch for the occupied states just below the gap 
edge at large ky^ and the peak in the energy dispersion of 
this state implies a van-Hove-like singularity contributing 
to the net M. 

Conlcusions. We showed that the magnetic proper¬ 
ties of singlet-triplet heterostructures are extraordinary 
sensitive to the type of triplet pairing, and to subdom¬ 
inant pairing channels. The dependence of the magne¬ 
tization on the phase distinguishes between helical and 
chiral pairing. Vice versa, our results allow control of the 
interface magnetization by the phase difference. We con¬ 


sidered an interface of high transparency and identical 
band structure on both sides. Relaxing this assumption 
will reduce the amplitude of the leaking superconducting 
component. In the first two cases we expect a reduction 
in the amplitude of M((/)), while in the last example the 
magnetization is weakly sensitive to the interface mis¬ 
match except near p ^ ir. Finally, our results suggest 
that such junctions act as possible spin pumps. A volt¬ 
age bias across the heterostructure would make the phase 
difference time-dependent, p = 2eVt, and the dynamics 
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FIG. 3. Spatial dependence of the spin-polarization. Panels 
(a)-(c) correspond to the cases of helical, chiral, and mixed 
symmetry order parameters as discussed in the text. Panel 
(d) shows a typical energy spectrum (for mixed symmetry 
case and 0 = 0). The thick line through the gap is the chiral 
edge state away from the interface. A significant contribution 
to the magnetization is due to the Andreev state just within 
the gap at large values of ky . 


of the magnetization is determined by this driving force 
as well as damping. If damping is small, for our last case 
of mixed-parity superconductor in contact with a con¬ 
ventional 5-wave system, the magnetization oscillations 
do not average to zero over a period of njeV. Detailed 
investigation of this effect is left for future work. 
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